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We make a short review about the synchronization in coupled phase oscillator models. Next, 
we study the common-noise-induced synchronization among active rotators. At an intermediate 
noise strength, the noise-induced synchronization takes place most effectively, which is analogous 
to the stochastic resonance. Finally, we study the synchronization of coupled phase oscillators with 
nonvariational interaction on scale-free networks. We find a sharp transition and a weak hysteresis 
in the nonvariational systems. The sharp transition is found also in the mean-field approximation. 
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Limit-cycle oscillation is a typical phenomenon in nonlinear-nonequilibrium systems. The synchronization or the 
frequency locking is a unique phenomenon in a population of coupled limit cycle oscillators^ Huygens observed 
mutual synchronization of pendulum clocks suspended from the same wooden beam in the 17th century. Buck and 
Buck reported synchronization of rhythmical flashes of fireflies in South Asia. 2 Recently, mutual synchronization of 
cellular activities was optically observed in the suprachiasmatic nuclei, which is the center of the circadian rhythm.— 
Rhythms such as the theta rhythms in the Hippocampus and the gamma rhythms in the visual cortex are considered 
to play some important roles in information processing in the brain4 

The theoretical study of the forced entrainment of a limit-cycle oscillator by an external sinusoidal force or the 
^Zh \ frequency locking between two limit-cycle oscillators has a long history in various research fields such as the classical 
mechanics and the electrical engineering. In a problem of circadian rhythms, Winfree first studied mutual entrainment 
in a population of coupled limit cycle oscillators. Kuramoto proposed a soluble coupled phase oscillator model, which 
. has a simple form and exhibits a phase transition analogous to the magnetic phase transition. 6 The Kuramoto model 
■ is written as 

i j r , N 

where fa and u>i represent the phase and the natural frequency of the zth oscillator, K is the coupling constant and 
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N is the total number of the oscillators. The natural frequency u>i is random, and the distribution of uj is denoted 
as g{oS), which is often assumed to be mirror-symmetric with respect to the average value ujq. Kuramoto showed 
that a phase transition from a disorder state to a mutually synchronized state occurs at a critical coupling strength 
K c . He showed that an order parameter a = (l/N) J^j exp(i</>j) for the phase order increases continuously from zero 

as a oc (A" — Kc) 1 / 2 , when the coupling strength K is slightly larger than K c . A macroscopic number of oscillators 
are entrained to the macroscopic oscillation, whose frequency is equal to the average frequency ujq of the natural 
frequencies. The ratio r of the entrained oscillators over the total number N also increases as r oc (K — Kc) 1 / 2 . The 
ratio r is interpreted as an order parameter for the frequency order. Strogatz and Mirollo discussed the stability of the 
disorder phase characterized by a = for K < K c - 7 Crawford derived an amplitude equation for the order parameter 
a(t) near k c & Recently, the Kuramoto model has been intensively studied and an review article of the model was 
written by Acebron et ali£ 

Various generalized coupled phase oscillator models have been studied. Sakaguchi and Kuramoto studied a nonva- 
riational model 
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where a is a parameter for the nonvariational coupling*^ This type of coupling characterized by nonzero a appears 
naturally in the general phase description of coupled limit-cycle oscillators ^ Wiesenfeld, Colet, and Strogatz derived 
this type model equation for coupled Josephson junctions?^ The coupled Joscphson junction model was further studied 
by Sakaguchi and Watanabe. 1 ? The model with a = is simpler and analogous to the XY model in magnetic systems, 
however, the nonvariational model with nonzero a is more natural as a nonequilibrium model. This generalized model 
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(2) is also soluble and exhibits a phase transition, however, the frequency of the macroscopic oscillation is generally 
larger than the average value u>q of the natural frequency. It is often observed in realistic systems that the actual 
frequency, which appears as a result of mutual entrainment, is larger than the average value of the natural frequencies. 
Daido studied even more generalized coupled phase oscillator models: 

^ = ^i + jj^Kfa ~<t>i), (3) 

and showed that the order parameter a can increase as a ~ {K — K c ), when h(ip) includes the second harmonics such 
as sirx(2 , i/;)»i* We proposed another soluble model including a noise termi^: 



5^sin(& -&)+&(*), (4) 



where is a Gaussian white noise satisfying (£,i(t)t;j(t r )} = 2D8i t jS(t — t'). The noise £i(t) for the ith oscillator is 
independent of the noise £j (t) for the jth oscillator in this model. We also studied the synchronization in a population 
of coupled active rotators: 

d<f>i K v—v 

— = Ui - bsiTMpi + — sax[4>j - (pi). (5) 

3=1 

When the parameter b is increased up to u>i, each oscillator behaves like a relaxation oscillator, which is characterized 
by the repetition of long-time slow motion and short-time fast motion. It is because the phase motion becomes slow 
near <f> = n/2. If the parameter b is further increased to be larger than w,, the phase locking occurs near = tt/2 and 
the oscillator changes into an excitable element. The relaxation oscillation and the excitability are also observed in 
various actual systems. In a population of these coupled active rotators, it was shown that the frequency entrainment 
takes place at several frequencies lu = ujq + nAwi. (n is an integer.) 15 
We further studied a slightly modified model equation: 

<M>- K N 

— ^ =uj- 6sin0i + F sin uj t + £i{t) + j^^sin^ - 0,-), (6) 
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in the parameter region u> ~ br^ In this model, external noises and a periodic force are further applied to the coupled 
active rotator model Eq. (5). When b ~ u>, each element is close to the excitable system. Excitable elements under 
a periodic force and external noises can exhibit the stochastic resonance, that is, the response to the periodic force 
is enhanced by the external noises of intermediate strength^ In our model, each element is further coupled with 
the other elements with equal strength. The stochastic resonance is expected to be enhanced by the cooperative 
interaction. Indeed, we found that the order parameter exhibits a limit-cycle oscillation with large amplitude at an 
intermediate strength of external noises, which is a phenomenon analogous to the stochastic resonance^ 

Global coupling is assumed in these coupled phase oscillator models, that is, each oscillator interacts with all the 
other oscillators with the same coupling strength. The global coupling makes a mean-field analysis possible, therefore, 
the models can become soluble in some cases. In realistic systems, such a uniform and global coupling is not plausible. 
Sakaguchi, Shinomoto and Kuramoto studied locally coupled phase oscillators on square lattices and cubic lattices and 
discussed the possibility of phase transitions in finite dimensional systems*^ Each oscillator interacts with oscillators 
on the nearest-neighbor sites as 

^ =Ui + K ^2 sin(<?!>j - (7) 
3'e-fCj 

where the summation is taken on the nearest-neighbor sites Ki of the ith site. We found that the phase order described 
by the order parameter <j does not appear even on the cubic lattices, but the frequency order denoted by the number 
ratio r of entrained oscillators seems to exist in the three dimensional system. Hong et al. studied numerically the 
finite dimensional model up to d = 6 using the finite-size scaling and found the lower critical dimension for the phase 
order is df = 4*^ 

Recently, complex networks such as the small-world network™ and the scale-free network^! have intensively been 
studied. In this trend, the coupled phase oscillators with variational interaction a — have been studied on complex 
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FIG. 1: (a) Average value of —bcos(j> at b = 1,Wq = 3 as a function of D 1 ^ 2 for the stationary distribution (12). (b) Temporal 
average of a for model (8) as a function of D 1 ^ 2 at K = 0.018 (circles) and K = (crosses) for N = 1000, to = 3 and 6=1. 
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FIG. 2: (a) Time evolution of a(t) for K = 0, u = 3, b = 1 and D 1/2 = 1.8. (b) Resultant frequency uj % = (cj>(t + T) - <j>(t))/T 
vs. toi with T = 1500. (c) Resultant frequency a), = {<j>{t + T) - 4>{t))/T vs. with T = 5000. 

networks such as small-world networks and scale-free networks. Hong et al. studied the synchronization on the 
small-world network^, and Moreno and Pacheco studied the synchronization on the scale-free network, and they 
found a desynchronization-synchronization transition similar to the one found in the globally coupled system^. 
Ichinomiya suggested that the critical coupling strength K c becomes zero in random scale-free networks with link 
number distribution P{1) oc Z~ 7 of 2 < 7 < 3^ Lee studied the system with the mean-field method^ Hong et al. 
studied finite-size scaling near the critical point on complex networks. 26 

In this paper, we study two new coupled phase oscillator models, which are closely related to the previous models. 
One model is globally coupled active rotators under common external noises, and the other is coupled phase oscillators 
with nonvariational interaction of nonzero a on a scale-free network. 



II. COMMON-NOISE-INDUCED SYNCHRONIZATION 



Various noise effects such as the stochastic resonance have been intensively studied*^ It is known that common 
noise can induce complete synchronization even for uncoupled homogeneous oscillators! 27 ' 28 Mori and Kai observed 
noise-induced entrainment in brain waves of humans^ In their experiment, mutual synchronization of a waves was 
induced by a common noise applied through one eye. We studied a one-dimensional active-rotator system with unstable 
coupling under a common noise, and found that mutual synchronization occurs by the common noise. The noise- 
induced synchronization was most facilitated at an intermediate noise strength, which is analogous to the stochastic 
resonance^, We found that desynchronization occurs intermittently near a critical parameter where the complete 
synchronization breaks down. In this section, we study a globally coupled active-rotator system with random natural 
frequencies under a common noise. The model equation is 

deb- K N 

— j = Wi - 6sin0 4 + — ^sin(0j - (8) 
where is a random natural frequency, the probability distribution of u>i is assumed to be a Gaussian distribution 
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with the average value loq and the variance k 2 , and the common noise satisfies = 2DS(t — t'). This 

model equation has a form of a mixture of Eq. (4) and (6), however, the common noise £(t) is applied to all active 
rotators, and therefore, the subscript i does not appear in the noise £(t) in Eq. (8). If two uncoupled active rotators 
with the same natural frequency are driven by the common noise, the equations of motion are written as 



~dt 

d<p 2 



- = cj - 6 sin 0i + £(i), 
at 



LU -bsm(f) 2 +t(t). (9) 

In the completely synchronized state, 4>i(t) = 4>2{t) is satisfied. The stability of the completely synchronized state is 
studied by a linearized equation: 

= -6cos0(i)<50, (10) 

dt 

where 6<f) denotes small difference: 5<fi — 4>2{t) — 4>i{t) and <j){t) denotes the completely synchronized motion: cf>(t) = 
4>i(t) = <f>2(t)- The long time average of —bcoscj) determines the linear stability for the synchronized motion. For 
D = 0, the time average of —b cos <fi is zero, and therefore the synchronized motion is neutrally stable. For nonzero D, 
it is possible that the long time average of —bcos(f> becomes negative and the synchronized motion becomes stable. To 
calculate the average of — bcoscj), we use the Fokker-Planck equation corresponding to the Langevin equation Eq. (9). 
For nonzero D, the probability distribution of 4> obeys the Fokker-Planck equation: 

dP d d 2 P 

- = -- { ( U0 -bBbxm+D W (id 

The stationary solution P s is given by2£ 

C ( e -2wu /D _i) r<t> e (-^oV>-f>cos «/>)/£> (W> 1 

P s {4>) = P s (0)exp{(^ -b + bco^)/D} 1 + ^ U ° — — , (12) 

I Jo e ^ J 

where P s (0) is determined by the normalization condition. The long time average of —b cos <f>(t) can be replaced by the 
ensemble average by the stationary distribution P s {4>). Figure 1 displays the average value of —bcos(f> as a function 
of D 1 / 2 for oj = 3 and b = 1. The average value of —6 cost/) is indeed negative, and it takes a minimum value at 
D 1 / 2 ~ 1.7, where the stability of the synchronized state is the strongest. In other word, the synchronization by the 
common noise is expected to takes place most strongly near D 1 ' 2 = 1.7. 

We have performed direct numerical simulation of the coupled phase oscillator model (8) with N = 1000, b— l,u>o = 
3 and k = 0.01. Figure 1(b) displays the temporal average of a — |l/iV 5Zj=i ex P(*^j)l as a function of D 1 / 2 for K = 
and K = 0.018. The order parameter represents a degree of the synchronization among the phase oscillators. The order 
parameter takes a maximum value near D 1 / 2 ~ 1.7, as is expected from the stability exponent. These simulations show 
that the common noises induce the synchronization among active rotators even if the natural frequencies arc randomly 
distributed. The synchronization is most effective at an intermediate strength of noise, which is also analogous to 
the stochastic resonance. The value of the order parameter takes a larger value for K = 0.018 than K — for any 
D, because the effect of the mutual entrainment is further added in case of nonzero K . However, the noise induced 
synchronization is more clearly seen in the uncoupled system K = 0, in that the difference of the maximum and the 
minimum values of the order parameter as a function of D is larger for K = than K = 0.018. Figure 2(a) displays 
a time evolution of a(t) for K — and D 1 / 2 = 1.8. The time evolution is strongly intermittent. It is because the 
stability exponent —bcoscp fluctuates in time markedly owing to the common noise. That is, — 6 cos </>(£) is negative 
on the average and the order parameter tends to take a nonzero value, however, it is possible that — b cos (f>(t) takes a 
positive value for a rather long duration, then, the synchronized motion becomes unstable and the order parameter 
is decreased to nearly zero. Figures 2(b) and (c) displays a relation of the numerically obtained resultant frequency 
Q l = (0(t + T)- ( /) l (t))/T ((b) T = 1500 and (c) T = 5000) vs. to, for K = and D 1 / 2 = 1.8. Step structures including 
many fine steps are seen in these figure, which imply mutual frequency entrainment at many different frequencies 
for a finite value of T. Thus, we have shown using this model that the noise-induced synchronization occurs even in 
a population of active rotators with random natural frequencies. The step structure becomes finer and finer as T is 
increased as shown in Figs. 2(b) and (c) and the structure disappears at T — > oo, because the complete frequency 
entrainment is impossible. However, the phase synchronization characterized by nonzero a(t) is maintained as shown 
in Fig. 2(a). 
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III. COUPLED PHASE OSCILLATORS WITH NONVARIATIONAL COUPLING ON A SCALE-FREE 

NETWORK 

In this section, we study a population of coupled phase oscillators with nonvariational coupling on a scale-free 
network. The model equation is written as 

^ = Ui +K y (sin(<£j - (j>i- a) + sin a), (13) 

where each oscillator is set on a scale-free network and the summation of the mutual interaction is taken only between 
oscillators connected by links in the scale-free network. We studied such nonvariational systems on linear, square and 
cubic lattices, and found that phase waves such as target waves and spiral waves appear inside of mutually entrained 
domains in the finite-dimensional nonvariational systems^! The nonvariational coupling is essential for the formation 
of the phase waves. The phase waves play very important roles for mutual entrainment. As a result, the mutual 
entrainment in systems with nonzero a can be rather different from that in the system of a — 0. For example, we 
found the multi-stability of mutual entrainment in a system of a — 7r/4 on a square latticed 

Moreno and Pacheco^ studied coupled phase oscillators with variational interaction of a = on a scale-free 
network proposed by Barabasi and Arbert^ The scale-free network by Barabasi and Arbert has a power-law degree 
distribution P(l) ~ l~ 7 with 7 = 3. Moreno and Pacheco found that the order parameter increases continuously from 
zero as K is increased. They analyzed the system by a finite-size scaling method and found that a phase transition 
occurs at a finite K. The order parameter increases a ~ (K — K c )@ with /3 = 0.46, which is consistent with the 
critical exponent 0.5 in the original globally coupled system (1). However, Lee suggested that the exponent (3 changes 
as (3 = 1/(7 — 3), which becomes 00 at 7 = 3. We have performed similar numerical simulations of coupled phase 
oscillators with nonvariational interaction for N = 2000. The preferential attachment algorithm by Barabasi and 
Arbert is used to construct a scale-free network. There is a free parameter m in the Barabasi- Albert model. The 
parameter m denotes the number of older sites to which a newly generated node attaches. In a scale-free network 
with parameter m, the average link number is expressed as (I) ~ 2m. We have first performed numerical simulations 
for a scale- free network with m = 2, and the probability distribution of uj is assumed to be a Gaussian distribution of 
k = 0.1 and luq = 0. 

Figure 3(a) displays a as a function of K for a — 0,7r/12,7r/6 and 7r/4. The coupling constant K is stepwise 
increased in the numerical simulation, that is, a final state of <j>i at a coupling constant K is used as a fresh initial 
condition at a slightly larger coupling constant K + AK. For a = 0, we found a continuous transition which is similar 
to that investigated by Moreno and Pacheco. However, the phase transition becomes sharp at a = 7r/12,7r/6 and 
a = 7r/4. The critical values are evaluated respectively at K c — 0.039,0.071 and 0.114 for a = 7r/12, 7r/6 and 7r/4. 
Figure 3(b) displays a as a function of K for a = 7r/6 and 7r/4, however, in this simulation, the coupling constant K 
is stepwise increased, and then K is inversely stepwise decreased. For these parameter values of a, the transitions are 
discontinuous and a hysteresis appears, although the bistable range is rather small. The bistability was not clearly 
observed for a = 7r/12. We have checked that the sharp transitions occur also for other scale-free networks with 
m = 2, 3, 5, and 8. Figure 3(c) displays the order parameter a when the coupling strength K is changed stepwise at 
a = 7r/4. (Here, a different scale-free network from the one used in Figs. 3(a) and (b) was constructed for m = 2.) 
We can evaluate the critical coupling strength K c as 0.1105,0.0625,0.033 and 0.0195 respectively for m = 2,3, 5 and 
8. As m is increased, (l)K c = 2mK c seems to approach a constant value about 0.29. 

In the scale-free network, there are several hub sites which might play an important role even for the mutual 
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FIG. 4: (a) Resultant frequency lo vs. link number U at m — 2, a — tt/6 and K — 0.07 in a desynchronized state, (b) Resultant 
frequency Q vs. natural frequency uj in the desynchronized state, (c) Resultant frequency uj vs. link number U in a synchronized 
state, (d) Resultant frequency uj vs. natural frequency uj in the synchronized state. 
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FIG. 5: (a) Order parameter a as a function of K numerically obtained from the self-consistent equation (16) at a — (dotted 
curve), a — tt/6 (dashed curve) and a = tt/4 (solid curve) for re = 0.1 and Pi(l)- (b) Order parameter a as a function of K at 
a — (dotted curve), a = n/6 (dashed curve) and a — 7r/4 (solid curve) for re = 0.1 and Pa(0- 



entrainment. In the coupled oscillators with nonvariational interaction a > 0, oscillators with larger connections tend 
to have larger frequency, if mutual synchronization does not occur. It is because the interaction term of the right-hand 
side of Eq. (13) is decomposed into 

K {sin((f>j — 4>i — a) + sin a} = K cos a sin(</>j — fc) + K sina{l — cos(<f)j — <fii)-} (14) 

The second term of the right-hand side of Eq. (14) makes the resultant frequency increase. Because oscillators at 
hub sites have a large number of connections and the resultant frequencies@increase owing to the summation term, 
which makes the mutual entrainment difficult at small K for the hub oscillators. However, for sufficiently large 
K, mutual synchronization proceeds and the phase differences \<pj — <f>i\ become small, and the effect by the term 
^sincv{l — cos(0j — 4>i)} becomes weak, then the mutual synchronization rapidly extends via the hub oscillators. It 
might be a reason of the bistability and the hysteresis. Figures 4(a) displays a relation of the link number U and the 
resultant frequency uii for each oscillator, and Figure 4(b) a relation of the natural frequency LOi and the resultant 
frequency u>i at a — tt/6 and K = 0.07 in a desynchronized state with a ~ 0.13. The scale-free network of m = 2 is 
used. It is clearly seen in Fig. 4(a) that the resultant frequencies are large for oscillators with large link number. No 
clear frequency locking is seen in Fig. 4(b). Figure 4(c) and (d) display a relation of li and a),-, and a relation of oji 
and bJi at the same parameter a = tt/6 and K = 0.07 in a synchronized state with a ~ 0.58. In this synchronized 
state, the oscillators with large link number are well entrained to the frequency Q ~ 0.0408. The resultant frequency 
uj ~ 0.0408 is definitely larger than the average value ujo = of the natural frequency. It is characteristic of the 
oscillator system with nonvariational interaction. There is a tendency that the phase <pi for the oscillators with larger 
link number I is larger in a population of the mutually entrained oscillators. That is, the hub oscillators tend to be 
in leading positions of the oscillation or play a role of pacemakers among the mutually-entrained oscillators. This is 
analogous to the fact that the oscillators with large natural frequencies play a role of pacemakers, which appear as 
centers of target waves, in finite-dimensional systems with nonvariational interactional 

The mean-field analysis might be applicable also in this system, if the phase correlation owing the network structure 
is neglected ! 10 i 24 ' 25 i 26 The order parameter weighed with the link number is defined as 
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where P(l) is the probability distribution of the link number I, another phase ip is introduced as ip = cp — Qt + a, and CI 
is the frequency of the order parameter. The phase ipi(l) of the ith oscillator with link number I obeys approximately 
the equation: 

— = u>i — fl + Kl sin a — Kla sin ip^. (15) 

If \u>i — lu + Kl sin a\ is smaller than Kl, the oscillator is entrained to the mean-field and the stationary solution of 
Eq. (15) is = sin - {(wj — uj + Kl sina)/Kl}. Since the natural frequency obeys the probability distribution 

g(oj), the phase distribution n s (ip; I) of the entrained oscillator is expressed as 

n s (ip; I) — g(fl — Kl sin a + Kla sinip)Kla cos ip, \ip\ < ir/2. 

If \u>i — uj + Kl sma\ is larger than Kl, the oscillator is not entrained to the mean-field, then the phase distribution 
ndsi^'j w, I) of the desynchronized oscillators is expressed as 

nd s (tp] uj,1) oc 1/\uj — f2 + Kl sin a — Kla sinip\. 

The order parameter can be expressed using the two phase distributions n s and rids- Then, we get a self-consistent 
equation for the order parameter and the frequency Q, as 

ae ia p ( 1 ) 1 = p ( l ) Kl2(7 jy 7 dipg{Q- Klsina + Klasiaip) costpe^ + iJ(l)^ , (16) 

where 

f 7 ^ 2 cosiifl — cos^~) 

J(0 = / # —{g{Q- Klsbxa + Kla/ simp) - g{VL - Kl sin a - Kla / siniP)} . 

Jo sin J ip 

If a = 0, O = and the critical value K c is calculated as K c — 2(1) /{Trg(0)(l 2 }}. If P(l) oc l/l 3 as in our scale-free 
network, K c become zero because (l 2 ) oc X)ti(^ 2 /^ 3 ) = 00 ■ However, if the system size is finite, the maximum number 
l m of the link number is finite. In that case, the critical value K c is not zero, and it is not so small because (l 2 ) 
diverges weakly as lnl m . If a is not zero and a is assumed to be infinitesimal, Eq. (16) is reduced to 

cosa(Z) = K^P{l)l 2 {n/2)g{Vl- Klsuva), 

sina(Z) = K^P{l)l 2 J^VL - Kl sin a), (17) 

where Ji(f2) = J °° <&e{<7(S! + x) — g(Q — x)} /(2x)*& The critical values K c and il c are obtained as a solution of the 
coupled equations (17). 

We have numerically solved the self-consistent equation (16) using two link number distributions, i.e., Pi (I) calcu- 
lated from the numerically constructed scale-free network for N = 2000 and m = 2 used in the previous numerical 
simulations shown in Fig. 3 and 4, and oc l/l 3 with l m = 1000. The natural frequency distribution g(uj) is 

assumed to be the Gaussian distribution with k = 0.1 and ojq — 0. The average link number (I) is 3.94 for Pi(l) and 
1.37 for i^CO- The self-consistent equation was solved by an iteration method. Figures 5(a) and (b) display a as a 
function of K at a = (dotted curve), 7r/6 (dashed curve) and a = 7r/4. (solid curve) for (a) and (b) -P2(0- 

For a = 0, continuous transitions from the desynchronized state to the synchronized state occur at K c = 0.0037 for 
Pi (I) and K c — 0.035 for P^Q)- The critical values and the order parameter a sensitively depend on the link number 
distribution. Here, the distribution P2G) decreases monotonously as l/l 3 . However, in the numerically constructed 
network characterized by Pi(l), the largest hub site has a link number U = 430 and the second-largest hub site has a 
link number U = 249, although Pi is roughly proportional l/l 3 for / < 20. For the exponent (3 in a oc (K — K c )° 

is rather large as was suggested by Lee^, but p3 is about 0.5 for Pi (I), which is consistent with our direct numerical 
simulation and the numerical simulation by Moreno and Pacheco. For a = tt/6, discontinuous transitions from the 
synchronized state to the desynchronized state occur at K = 0.047 for Pi(l) and at K = 0.128 for P2{1), when K is 
continuously decreased. For a = 7r/4, discontinuous transitions from the synchronized state to the desynchronized 
state occur at K = 0.065 for Pi(l) and at K = 0.128 for P2(l), when K is continuously decreased. Transitions from 
the desynchronized state to the synchronized state occur respectively at K = 0.059 and K = 0.079 for a = tt/6 and 
7r/4 in Fig. 5(a), and at K = 0.134 and K = 0.149 for a — tt/6 and 7r/4 in Fig. 5(b), when K is continuously increased. 
That is, hysteresis is seen also in the mean-field analysis. For a = tt/A and Pi (I), the average of the two critical values 
is estimated as K c ~ 0.072, then, (l)K c ~ 0.285, which is close to the critical value for large m in Fig. 3(c). These 
results of the mean- field approximation for Pi(l) are qualitatively consistent with the results of the direct numerical 
simulation shown in Fig. 3. However, quantitative agreement is not so satisfactory, probably because the mean-field 
approximation is not exact on the scale-free network and the finite-size effect is rather strong. 
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IV. SUMMARY 



In this report, we have made a short review about the Kuramoto model and the generalization. Next, we have 
studied the common-noise-induced synchronization among active rotators. At an intermediate noise strength, the 
noise-induced synchronization takes place most effectively, which is analogous to the stochastic resonance. Finally, we 
have studied the synchronization of coupled phase oscillators with nonvariational interaction on scale-free networks. 
We have found sharp transitions and weak hysteresis in the nonvariational systems. The mutual entrainment via 
hub oscillators seems to be important in these nonvariational systems. We have numerically solved the self-consistent 
equation by the mean-field approximation. We have found that the discontinuous transition occurs also in the mean- 
field approximation for nonzero a. 
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